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Abstract 

^ ! A chaotic model of the early universe within the framework of the singularity- free 

a^ : 

^ ■ solutions of Einstein's equation is suggested. The evolution of our universe at its 

» ' 
Q . 

Q . early stage, starting out as a small domain of the parent universe, is governed 

O ■ 

CN I by the dynamics of a classical scalar field (j) . If in any such domain, larger than 

^ ■ Planck length, happens to be very large, (/> may develop a dominant inhomogeneous 
00 ; 

^ \ mode, leading to an anisotropic inflation of the universe. The particle 4> is coupled to 
O ; 

^ ■ other particles, which are produced copiously after inflation and these thermalize 

O . leading to a rather low temperature universe (T > 10^ Gev). The electroweak 



B-(-L Baryogenesis is assumed to account for the observed baryon asymmetry. The 
universe now passes through a radiation-dominated phase, leading eventually to a 



X 

H 

5r ! matter-dominated universe, which is isotropic and homogeneous. The model does 
not depend on the details of Planck scale physics. 
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In the chaotic inflationary models, the universe is assumed to evolve out of a random 
distribution of initial data. This makes the model superior to other inflationary 
models, which require a very special class of potentials and/or speciflc initial 
conditions. For example, in the new inflationary model, the universe is assumed to 
be in thermodynamical equilibrium initially, with the inflaton field (/> at the local 
minimum of its temperature - dependent potential V(0,T). On the other hand, the 
chaotic model should work, in principle, for a random distribution of the initial 
values of (j) and for a large class of potentials. However, as Linde ^ has pointed out, 
the values of (f) are not totally unspecified. When one specifies the initial conditions 
for a classical space-time, one must have 

p < V{4>) < {dicj>f < M^, i = 0, .., 3 (1) 

These constraints make the scalar field (f) almost homogeneous, if it is large, as 
can be seen from the following considerations. A closed mini-universe, emerging 
out of the quantum space-time foam, has a typical size of Planck length, M~^. 
In such a universe, with the gradient being less than Mp,the field (f) cannot differ 
from some average value 0o by more than Mp. Thus, if the value of 0o » Mp, (j) is 
essentially homogeneous to begin with. Linde has assumed a similar homogeneous 
and large value for even for an open universe. However, the fact that ^ Mp 
leads to a very low and 'unnatural' value for the coupling constants (e.g. A ~ lO"-*^^ 
for a A(/)^ theory). Moreover, the situation with the open universe is not very 
clear. It would indeed be useful if a scenario with 0o < Mp also works. Obviously, 
one has to encounter an anisotropic and inhomogeneous space-time and also an 
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inhomogeneous </>. But such a model will be more 'natural' and also in keeping with 
the expectation that the universe would have evolved to its present state even if 
the initial geometry is not carefully chosen. A first step in this direction will be to 
consider a mini-universe with lesser symmetry (say one with cylindrical symmetry) 
and try to build a chaotic model with initial values of (j) less than Mp. We present 
here such a model, which differs considerably from the model of Linde. Also, in 
our model, the universe need not pass through states of very high temperature or 
density. The universe remains inhomogeneous and anisotropic till very late (till 
matter- dominated era). The inflation takes place in two stages. Our description of 
space-time makes use of the singularity- free solutions^ These solutions describe 
an inhomogeneous and anisotropic space- time, containing a perfect fluid. These 
may be generalized to describe different stages of the evolution. Ruiz and Senovilla^ 
have given a useful factorizable expression for the metric, describing a space-time, 
which admits two commuting space-like Killing vector fields, which are mutually 
as well as hypersurface orthogonal. We follow their notations but consider only a 
subfamily, viz. 



ds^ ^T''+^ F^{-dt'^ +H^dr^) + T^+^'GPde^ + T^-'^GP-^dz^ (2) 



where T = T{t) satisfies the equation 



^ = ea^e=0,±l, (3) 



with a an integration constant. The functions F, H, G, P are functions of r only. 
These are related by 



= (P/G)" (4) 
and H, G and P satisfy the differential equations : 

H' 

a' + aj3 — a — = ena^H^ (5) 
H 

TTl 

(1 - n)(3' + - 2n(3a + - (1 - n) /? — = n(l - n)ea^H^ (6) 

H 

with q; = -p and j3 — ^ . In above, dot and prime denote differentiation w.r.t. 
t and r respectively. Since there are only three relations among the four functions 
F,H,G,P, they will be given in terms of an arbitrary function C(r). However, we shall 
impose, when convenient, an additional relation, thus determining these functions 
completely. 

Some interesting features of these non-singular solutions have already been 
studied^-^ The singularity free family (2) can be considered as arising out of FRW 
metric for open universe through a natural inhomogenisation and isotropisation 
process^ and it is unique^ and geodesically complete^. If one wishes to study the 
cosmology of a closed universe, evolving out of the quantum foam, the non-singular 
nature of the solution is of no special advantage. However, we wish to study open 
universes as well. These non-singular solutions may be useful, if one tries to give 
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a classical description of an ever-existing parent universe. We will consider this 
problem elsewhere. For the purpose of the model, it is sufficient to assume that there 
is an ever-existing parent universe and to follow the evolution of a small domain 
(larger than Planck size) of this universe, which happens to satisfy suitable initial 
conditions, to be given later on. The evolution of this mini-universe is governed by 
the dynamics of a scalar field (/>, having an effective potential V(</)). It can be shown 
that the assumed cylindrical symmetry of the space-time ensures that cf) will not 
depend on 6 and z, i.e.0 = r). The relevant field equation, with the metric (2), 
is given by 

Note that the equation gets simplified, if we impose the condition G = H. We shall 
follow the evolution in steps and outline the relevant features : 

1. First inflationary stage : 

In this model, in order to qualify as a mini-universe, a small domain of space-time 
has to satisfy the following conditions : 

Initially , should be much larger than as well as H = G), so that the 

relevant equation of state is that of stiff matter. 

When these conditions are met, the mini- universe can have a metric of the form 
(2), for n < — 1. With G = H, we have a class of solutions 



C(r) = [-na^M{r + q)]-'', 
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P^(r) = n^a^M^C^{r), H{r) = G{r) = 



a{r + q) ' 

T{t) = ToCosh{at), XP-XP= ^(/>'T-("+i) C''^ (8) 

where M, q and a are constants. We note that the above solution leads to an 
inhomogeneous expansion given by 

V{t) = V{0)exp[- ^^"l^J ■ T(t)-^("+i)] (9) 
{n+ 1)1' [r) 

where V(0) is the volume of the miniuniverse at t = 0. The solution with n = — 3 
may be useful to describe an ever-existing classical space-time or a part thereof, as 
it permits no expansion. While any solution with —3 > n > — 1 is relevant, we 
consider, in particular, the solution with n = — 2. It gives an inflationary universe 
which expands as 
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V{t) = V{0)exp[-—^cosh{at)]. (10) 
r [r) 

The function F(r) in this case is given by 

F(r) = (11) 

so that a suitable choice for q and a can provide adequate inflation everywhere. The 
space-time is highly anisotropic, there being a rapid expansion along the z-axis and 
contractions along directions normal to it. The inflation ceases because the fleld (j) 
increases rapidly. To see this, we rewrite equation (7), neglecting the (f)" term, as 
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^ + 4f^ = -4A*^T'*+\ (12) 

where we have substituted V(0) = and ^ = (j)F{r). Since is small initially, 
we neglect the R.H.S. and obtain a solution for the increasing mode, 

* = *o + — ton-He'^O - ^ (13) 

where is the value at t = 0. We choose A as a constant and not a function 
of r to restrict the discussion to the case of simple inhomogeneity. Note that ^ 
approaches rapidly to a limiting value, \E'i = \E'o + ^ and if \E'o everywhere in 
the miniuniverse, we have almost a homogeneous ^ . However, this should happen 
before the potential term becomes important and the stiff matter equation of state 
gets appreciabely disturbed. This gives the condition 



or 



2a^ 

We, however, want ^ to be much smaller than Mp. The singularity free class 
of solutions seem to be very rich and can possibly accommodate the case when 
the potential term makes a contribution to the energy density and pressure. The 
class also includes a solution relevant for a cosmological constant {m — — —1) 
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but the inhomogeneous model does not pass through this stage. Thus, the usual 
potential-driven inflation does not occur in this model. As the held increases, 
there is a change in scenario, because is presumed to be unstable, being coupled 
to other flelds, both scalars and fermions. Eventually, (f) will decay producing 
these particles, which thermalize. It is known^ that one can get a hierarchy of 
temperatures depending on how </> decays. The mechanism of baryogenesis gives 
a lower bound on the temperature. If baryogenesis is due to B + L violating 
electroweak processes, a ' reheating' temperature ~ 10^ Gev is adequate. There 
does not seem to be any compelling reason for the universe to pass through a stage 
of higher temperatures. After thermalization, the universe is full of radiation and 
another inflationary stage begins. 

2. Second Inflationary stage : 

A radiation dominated universe can be described by the metric (2) for n = 3. This 
is the original Senovilla solution. We have indeed a large class of solutions : 

ds^ = T^C^i-dt'' + H^dr^) + T^CPde"^ + T-^GP'^dz^ (15) 
with T = Tocosh{at),C{r) an arbitrary function, 

p = 3p=^k T-^C-^ (16) 

with k ^ 0. There is now a contraction along the z- direction and expansions along 
directions normal to it. This peculiar two stage inflation helps to make the universe 
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isotropic. The expansion during this stage can be calculated from the relation, 



with t2 > ti. Since both energy density and pressure fall as , it is obvious that 
massive particles will soon become non-relativistic and a matter-dominated stage 
will be ushered in. 

3. Matter - dominated era : 

After the second inflation, we expect the universe to become isotropic and 
homogeneous and go over to the FRW universe. Actually, within the framework of 
the metric (2), there is only one window through which this can happen and this is 
during the matter-dominated era. The solution is 




(17) 



n = 0, F(r) = H{r) = 1, G{r) = P{r) = r 



(18) 



corresponding to the metric 



ds^ = T{t)[-dt^ + cZr^ + r'^de'^ + dz% 



(19) 



Einstein's equations now reduce to 



3 

G'oo = T-i.-^=XP 



(20) 
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Gil = G22 = Gss = T-'[-^ + = 



(21) 



which has the solution, 



T{t) = Bt\p{t) 



12 



(22) 



Bt^' 



We can rewrite the results in the conventional form by introducing a new time rj , 
viz. 



which is the well-known Einstein-deSitter universe. The model thus leads to a flat 
space with p ~ a~^, a ~ r/^ . It explains in a natural way why the universe is so flat 
or why p is so close to the FRW critical density. We have introduced a new time rj 
, but this is the only time of physical interest, because no astrophysical phenomena 
of earlier epochs are observable to us. 

It is useful to show that the universe does tend to roll down to the symmetric state 
given by equation (19). A simple way to do this is to study the stability of the 
solution (19). We put H = 1, which can be done without any loss of generality. If 
we now consider small linear perturbations of the functions P(r) and G(r) about the 
solution (19), making use of the fleld equations (5) and (6), the solution is seen to 
remain invariant, modulo a trivial scaling. True to our expectation, the symmetric 
state is most stable, and given sufficient time, the universe will roll into it. 



ds^ = -dri^ + a'^{ri)[dr^ + r^dO^ + dz^] 



(23) 
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To summarize, we have presented a scenario involving a set of initial data 
not considered in Linde's chaotic model. In a realistic chaotic model of the 
early universe, one should consider the entire phase-space and also include the 
corrections due to quantum fluctuations of the scalar field. Although our model 
is based on a particular class of inhomogeneous solutions, it is quite possible that 
initial inhomogeneities of many different types may also permit evolutions leading 
eventually to a large homogeneous and isotropic universe. An understanding of all 
these possibilities will be useful to arrive at the Grand chaotic model. The effect of 
quantum fluctuations of cj) is presently under study and will be reported elsewhere. 
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